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Abstract 

The response of the Neuberger-Dirac fermion operator D = I + V 
in the topologically nontrivial background gauge field depends on the 
negative mass parameter 777.0 in the Wilson-Dirac fermion operator Dw 
which enters D through the unitary operator V = D W {D\ 1] D W )^ 1 ^ 2 . 
We classify the topological phases of D by comparing its index to the 
topological charge of the smooth background gauge field. An exact 
discrete symmetry in the topological phase diagram is proved for any 
gauge configurations. A formula for the index of D in each topological 
phase is derived by obtaining the total chiral charge of the zero modes 
in the exact solution of the free fermion propagator. 



PACS numbers: 11.15.Ha, 11.30.Fs, 11.30.Rd 



1 Introduction 



It is well known that the Neuberger-Dirac operator [[[]. § 

D = l + V, V = D^DlDv,)- 1 ' 2 (1) 

reproduces the exact zero modes with definite chirality and realizes the Atiyah- 
Singer index theorem on a finite lattice for smooth background gauge fields 
with nonzero topological charge || . However, this is only true for the negative 
mass parameter m approximately in the range (0, 2r w ) for the Wilson-Dirac 
fermion operator D w with Wilson parameter r w > 

D v , = -m + -[^{V*+V li )-r v ,V*V ll ] (2) 

where and V* are the forward and backward difference operators defined 
in the following, 

V^{x) = Un(x)ip(x + ft) -ip(x) 
V*ip(x) = ip(x) - Ul(x - p)ip(x - p) 

Except for the region m < which corresponds to the positive mass in D w 
that gives zero index of D for any background gauge field has been noted by 
Neuberger in ref. ||, there is literally no discussions in the literature about 
what would happen to D for other values of mo, namely, mo > 2r w . In ref. 
P]], it is remarked that D w must have an appropriately tuned negative mass 
term such that D is in the right phase to reproduce QCD. However, there is 
no further discourse on how to identify the phase such that D can possess 
the proper topological characteristics and reproduce the correct continuum 
physics. 

In this paper we investigate the topological phases of D with respect to the 
negative mass parameter m and the background gauge field respectively, by 
comparing its index(D) to the topological charge Q of the background gauge 
field. If index(D) = for any smooth gauge configurations, then D is called 
to be in the topologically trivial phase. If index(D) = Q, then D is called 
to be in the topologically proper phase. If D is not topologically trivial but 
index(D) ^ Q, then D is called to be in the topologically improper phase. 

First, we sketch schematically the topological phases of D in Fig. [|. On 
a d- dimensional lattice^], the topological phases of D in the free fermion limit 
divide mo values into d + 2 intervals according to the real eigenmodes of D. As 
we turn on the background gauge field, the locations of the phase boundaries 
also evolve accordingly. The first region cuts the interval mo < in the free 
fermion limit. In this phase, the spectrum of D does not possess any zero 

1 In this paper we always assume that d is a positive even integer. 
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modes and thus cannot describe massless particles, and the index is zero for 
any background gauge field. Therefore D is called topologically trivial in this 
phase. The second region cuts the interval < m < 2r w in the free fermion 
limit. In this phase, D has exact zero modes with definite chirality, and the 
index of D is equal to the topological charge Q top of the background gauge field. 
Therefore D is called topologically proper in this phase. The third region cuts 
the interval 2r w < mo < 4r w in the free field limit. In this phase, D also has 
exact zero modes with definite chirality, but the index of D is (1 — d)Qt op , not 
equal to the topological charge of the background gauge field. Thus D is called 
topologically improper in this phase. The next phase again has width 2r w in 
the free fermion limit. The process of dividing mo into different topological 
phases continues until it reaches mo = 2dr w in the free fermion limit, then the 
last phase is the region with the edge, 2dr w < m < oo. In any dimensions, 
only the region cutting the interval < m < 2r w in the free fermion limit 
is the topologically proper phase for D, other phases are either topologically 
trivial or improper. The index along each one of the phase boundaries starting 
at ( m = 2kr w , k = 0, • • • , d ) is zero identically for any background gauge 
field. 

In the free fermion limit, we can distinguish different phases of D according 
to its real eigenmodes. The phase transition is signaled by the change of the 
numbers of real eigenmodes. We note that in the free fermion limit the real 
eigenmodes of D must be either at the origin or the corners of the Brillouin 
zone. However, the converse of this statement is not always true, since at some 
values of mo, some of the eigenmodes at the corners of Brillouin zone become 
complex. Then a formula for the index of D in each phase can be derived ( in 
Section 3 ) by obtaining the total chiral charge of the zero modes in the exact 
solution of the free fermion propagator. For smooth background gauge fields, 
this formula ([34]) gives the correct index in each phase except at the phase 
boundaries ( m = 2kr w , k = 1, • • • , d ). 

For any arbitrary gauge configuration, we prove that there exists an ex- 
act reflection symmetry in the topological phase diagram with respect to the 
symmetry axis mo = dr w , i.e., index[D(mo)] = —index [D(2dr w — mo)], for 
a finite lattice of even number of sites in each dimension. It follows that if 
there is a phase boundary located at m = x < dr w on the LHS, then there is 
another phase boundary located at m = 2dr w — x on the RHS, and these two 
phase boundaries appear as images to each other with respect to the mirror at 
mo = dr w . Since d is an even integer, the exact reflection symmetry implies 
that the symmetry axis m = dr w is a phase boundary with zero index and its 
location is invariant for any background gauge configurations. ( for example, 
see Fig. 2 and Fig. 3 ) 

As we turn on the background gauge field smoothly from zero to some finite 
values, the locations of the phase boundaries also evolve accordingly. When the 
gauge fields reach a certain level of magnitude and roughness, the phase bound- 
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aries bifurcate and the phase diagram becomes very complicated. We expect 
that for gauge configurations in dynamical simulations the phase boundaries 
may change significantly from those of smooth background gauge fields. How- 
ever, the exact reflection symmetry must hold in general and provides a useful 
scheme to identify the topological phases especially for rough gauge configu- 
rations. We note that at one of the topological phase boundaries of D, the 
theory is free of species doublings in the free fermion limit but turns out to 
be topologically trivial in any background gauge fields. This provides an ex- 
plicit example that for D free of species doubling in the free fermion limit is 
not sufficient to gaurantee the realization of Atiyah-Singer index theorem on 
the lattice, in contrast to the condition of reproducing correct axial anomaly 
in perturbation calculations ||. This may indicate that the topological char- 
acteristics of the fermion operator D are of nonperturbative origins and thus 
could not be revealed by any perturbation calculations. We note in passing 
that the locality of Ginsparg- Wilson fermion operator D is not relevant to its 
zero modes and index, as discussed and demonstrated in ref. || [JJ. 
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Figure 1: The topological phase diagram of D on a (^-dimensional lattice 
for smooth background gauge field. The phase boundaries denoted by dashed 
lines are sketched to indicate their evolutions with respect to the background 
gauge field. The critical points in the free fermion limit ( at the horizontal 
line ) are determined exactly. The reflection symmetry with respect to the 
symmetry axis mo = dr w is also shown explicitly. 



The outline of this paper is as follows. In Section 2, we prove the exact 
reflection symmetry in the topological phase diagram of D. In Section 3, we 
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obtain the exact solutions of the eigenvalues and the fermion propagator of D 
in the free fermion limit, then use them to analyze the topological phases of 
D. A formula for the index of D in each phase is derived. In Section 4, we 
perform numerical experiments to investiagate the topological phase diagrams 
of D in two dimensional and four dimensional lattices respectively. In Section 
5, we summarize, and also demonstrate the bifurcation of topological phases by 
increasing the amplitudes of the local sinusoidal fluctuations of the background 
gauge field on a two dimensional lattice. In the Appendix, we discuss the 
spectral flow of the Neuberger operator as a function of the negative mass 
parameter mo, in particular, its behaviors at the phase boundaries, in order to 
gain a perspective on the underlying mechanism which gives the zero index at 
the phase boundaries. 



2 The discrete symmetry of the index 

In this section, we prove that for a finite lattice with even number of sites in 
each dimension ( L M = N^a, = even integer, fj, — 1, ■ ■ ■ , d ), the topological 
phase diagram of the Neuberger-Dirac operator D(mo) has an exact reflection 
symmetry with respect to the axis m = dr w , and the index of D is odd under 
the reflection, i.e., index[_D(m )] = — index [D{2dr w — m )}. 
Proof : 

For a given gauge configuration, the index of D(m ) is defined to be n_(m ) — 
n + (mo), where n^(n + ) denotes the number of zero modes of negative ( positive 
) chirality. Let 01, s = 1, • • • , n + be the normalized eigenf unctions of the zero 
modes with positive chirality, 

D(m )4> s + = 0, 750 s + = S + . (3) 

Then <p s + is a real eigenmode of D w with negative real eigenvalue — fi s such 
that 

Virno)^ = -0% (4) 
D w (m )(j) s + = -n s (j) s + (5) 

Similar equations for zero modes with negative chirality are understood. The 
standard Wilson-Dirac lattice fermion operator D w , Eq. (|2|) can be rewritten 
as 

[D w (m )] X! y = (dr w - m )5 Xty - [T(r w )] Xjy (6) 

where 

[T(r w )] x>y = - X] [(r w - l^U^x)5 x+flty + {r w + i^U^x - A)4- A J (7) 

Z A 4 
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The Dirac and color indices are suppressed. The fundamental properties of T 
are § : 

75T75 = rt (8) 

STS = -T (9) 

where S XiV = Ylu{~ l) x>1 $x,y which satisfies S 2 = 1 and = S. The first 
property is the usual hermiticity condition. The second property is also obvious 
except at the boundaries of the lattice ( with periodic boundary conditions ), 
where the forward and the backward difference operators in T would produce 
a factor +1 ( -1 ) for odd ( even ) number of sites in this direction. Therefore 
Eq. (|9|) is satisfied exactly on a finite lattice having even number of sites 
in each dimension, but only satisfied approximately on a finite lattice having 
odd number of sites in some dimensions. Of course, it is satisfied exactly 
on the infinite lattice. It is evident that an eigenfunction of D w must be an 
eigenfunction of T, and vice versa. Then from Eqs. (|5|) and (^|), we have 

T<j>% = t s <f) s + (10) 

where t s = // s + (dr w — m ). Using Eq. ([5J) and S 2 = 1, we obtain 

det(T-tJ) = det(S) det(T - tl) det(S) 

= det(STS - tl) = det(-T - tl) 

Thus the eigenvalues of T must come in ±t s pairs. It follows that the eigenval- 
ues of D w {mo) must come in dr w — m ± t s pairs. Now if we substitute m by 
2dr w — m , then the eigenvalues of D w (2dr w — m ) come in —(dr w — m =F t s ) 
pairs. Therefore, from Eq. (|5|), if _D w (m ) has a negative real eigenvalue — ji s 
for the eigenfunction <p s + , then D w {2dr w — m ) has a corresponding positive 
real eigenvalue, +fi s with the same eigenfunction (p s + , and this implies that 
D(2dr w — m ) has a +2 real eigenmode while D(mo) has a zero mode, and 
vice versa, i.e., 



nt{m ) 


= n + (2dr w 


- mo) 


(11) 


n 2 (m ) 


= n^{2dr w 


- m ) 


(12) 


n + (m ) 


= n,2(2dr w 


- m ) 


(13) 


n-(m ) 


= n2~(2dr w 


- m ) 


(14) 



As it was first proved in ref. || that the number of zero modes and the number 
of +2 modes for D{mo) must obey the relation 

n^(m ) - rq (m ) = n_(m ) - n + (m ), (15) 

we immediately obtain 

n + (2dr w — m ) — n^{2dr w — m ) = n_(m ) — n + (m ). (16) 
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That is 

index [D{2dr w — m )] = — index [D(m )]. (17) 

This is the discrete symmetry of index(D) which holds for any gauge config- 
uration. It transcribes to the reflection symmetry in the topological phase 
diagram with the symmetry axis at mo = dr w . When the gauge fields are very 
rough, the topological phase diagram becomes very complicated, the only thing 
surely survives in such chaotic environments is the exact discrete symmetry, 
Eq. (|ITD . For a finite lattice having odd number of sites in some directions, the 
reflection symmetry Eq. ( |17D is only an approximate discrete symmetry which 
tends to the exact reflection symmetry in the infinite lattice limit. One imme- 
diate consequence of the exact discrete symmetry ( |17D is that at the symmetry 
axis m = dr w , the index is zero for any background gauge field, 

mdex[D(dr w )\ = (18) 

Since d is an even integer, the symmetry axis mo = dr w is a phase boundary with 
zero index and its location is invariant for any background gauge configurations. 
( for example, see Fig. 2 and Fig. 3 ). 

We also note that T(r w ) has another discrete symmetry, 

[T(r w )]t = -T{-r w ) (19) 

which follows directly from Eq. ([?]). Then we have the following discrete 
symmetry for D w , 

[D w (rn , r w )] f = -D w (-m Q , -r w ) (20) 

From this relation, we immediately know that the zero modes of D(rriQ,r w ) 
are the +2 eigenmodes of D(—m , —r w ), and vice versa. Using Eq. fll5]) , we 
obtain the following discrete symmetry for the index of D 

index [D (m , r w )] = — index[D(— mo, — r w )]. (21) 

This discrete symmetry implies that if we change the sign of the Wilson pa- 
rameter r w , then the new topological phase diagram is the negative image of 
the old topological phase diagram with respect to the mirror at m = 0. This 
reflection symmetry will not be exploited in this paper. 

We remark that the discrete symmetries, Eqs. ( |17D and ([21] ) discussed 
above also hold for the generalized Neuberger-Dirac operator || 

D = 2M(1 + V)[(l -V) + 2MR(1 + V)]' 1 (22) 

where R is a non-singular hermitian local operator which is trivial in the Dirac 
space, M is an arbitrary mass scale, and V(m , r w ) is the unitary operator 
defined in Eq. ([!]). The reason for these discrete symmetries holding for the 
Ginsparg- Wilson generalized D in ( [Z2"D is that the zero modes of D are invariant 
for any R || [/J, and hence any relationships involving only the indices must 
be R- invariant. The proof for the general case is obvious from the above proof. 



6 



3 The exact solutions in the free fermion limit 

In the free fermion limit, the index of D must be zero, but we can use the num- 
ber of zero modes of D to distinguish different topological phases of D. Since 
there is a symmetry between the number of zero modes and the number of +2 
modes as first shown in ref. ||, Eqs. (pTTD -(|T4|) , then it suffices to use both the 
number of zero modes and the number of +2 eigenmodes as the order param- 
eter for the topological phases of D in the free fermion limit. Then a phase 
transition is signaled by any changes of these two numbers. Alternatively, we 
can examine the fermion propagator and use the number of massless ( primary 
and doubled ) fermion modes it contains to distinguish different topological 
phases of D. From the total chiral charge Q§ of the massless fermion modes 
in the fermion propagator, we can determine the integer multiplier of the FF 
term in the axial anomaly q(x;A,D) = tr[j5(RD)(x,x)] for D in a smooth 
background gauge field of topological charge Qt op , which also can be deter- 
mined by explicit perturbation calculations. Furthermore we can also predict 
that the index of D is equal to the total chiral charge times the topological 
charge of the background gauge field, i.e., 



This expression is deduced straightforwardly from the weak coupling pertur- 
bation theory, which in fact correctly predicts the index of D in smooth back- 
ground gauge fields for all values of m except at the phase boundaries starting 
at m = 2r w , 4r w , • • • , 2dr w , where the index is always zero regardless of the 
total chiral charge of the zero modes. The disagreement between the predic- 
tion ( |23| ) and the zero index at the phase boundaries, indicates that there may 
be some nonperturbative or topological effects which cannot be taken into ac- 
count in the perturbation theory. We will show that Eq. ( PB"D violates the 
exact reflection symmetry at the phase boundaries. 

First, we obtain the exact solutions of the eigenvalues and the fermion 
propagator of D. In the free fermion limit, the eigenvalues of the Neuberger- 
Dirac operator D [ Eq. (ffl) ] on a <i-dimensional lattice [ L M = N^a, \i = 
<i(even integer) ] with periodic boundary conditions are 



index(D) = Q 5 Q top 



(23) 




(24) 



where 



u{p) 



d 

m + r w ^2[l- cos^a)] 



(25) 




(26) 
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zzK^.^eC^e ifiV^iseven 

and each real X s has degeneracy 2^ consisting of chirality ±1 pairs in the 
Dirac space, while each complex X s has degeneracy 2[ _ 2 _ 1. 
The free fermion propagator in momentum space is |3j 

,(0)/ a - a 



2 «7 M sm(p M a) 



where 



T(p) = ±[tf(p)-u(p)] (29) 

The constant term a/2 which vanishes in the continuum limit can be ignored 
since its role is the same as the chirality breaking operator R in the Ginsparg- 
Wilson fermion propagator. The second term is equal to the naive lattice 
fermion propagator times the factor T(p). The suppression of the doubled 
modes is due to the last factor T(p) which becomes zero for some or all of the 
doubled modes depending on the value of mo- We have deliberately written the 
fermion propagator in the form of Eq. ( pip such that the decoupling mechanism 
for the doubled modes is the most appealing to us. It is interesting to note that 
the way doubled modes are decoupled in ( p8|) is very different from that of the 
Wilson fermion propagator. The former is a complete suppression even at finite 
lattice spacing, e.g., T(p) = for all doubled modes for < m < 2r w , while 
the latter is of order a which is completely decoupled only in the continuum 
limit ( a — > ). This implies that, the Neuberger-Dirac fermion can reproduce 
the correct continuum axial anomaly and the exact Atiyah-Singer index theorem 
at finite lattice spacing, while the Wilson fermion can only recover them in the 
continuum limit. In other words, any ( primary or doubled ) massless fermion 
mode in the Neuberger-Dirac fermion propagator emerges in a very clean-cut 
way, i.e., either existing or vanishing, for any finite lattice spacing. 

In order to discuss the chiral charge of the doubled modes, we consider 
the free fermion propagator in the continuum limit ( a — > ) on the infinite 
lattice with the momentum space, the Brillouin zone ( BZ ) which is specified 
by p^ G (—71 /a, n/a] for \i = 1, • • • , d. In the vicinity of the origin p = 0, we 
can expand sin(p M a) = p^a + 0(a 2 p 2 ^) and T(p) = m + 0(ap 2 ) in Eq. (f2~8|) to 
obtain 

SP(p) a ^°-^m + O(a 2 p) (30) 
p A 

where the constant term a/2 has been subtracted and the factor mo can be 
absorbed by using the generalized Dirac operator ( f22f ) with M = mo and R = 
l/2mo, thus D = mo(l + V). Then Sp\p) in (j30|) agrees with the free fermion 
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propagator in continuum. This is the primary mode of the lattice fermion. 
Besides the primary mode, it is also possible for Sp (p) to develop massless 
doubled modes at other regions of the BZ due to the presence of sinQo^a) rather 
than p^a in Eq. (|28| ) . It is evident that the massless doubled modes can only 
occur around the 2 d — 1 corners of the BZ, i.e., ®^ =1 {0, 7r/a} — {p = 0}. If the 
factor T(p) vanishes smoothly at any one of these corners, then the doubled 
mode at this corner is decoupled from the theory, otherwise this doubled mode 
will contribute to the propagator. Furthermore, a doubled mode may possess 
opposite chiral charge of the primary mode due to the property of sin^a). 
When p^ is near n/a, we must shift the momentum p^ to p^ = p^ — n/a 
such that the higher order terms 0(p 2 a 2 ) in sin(p^a) can be neglected in the 
continuum limit a — > 0. Then the relationship sinfj^a) = — sin(p a) produces 
an extra minus sign which can be absorbed by redefining 7^ = —7^. Since 
75 = ria=i 7a» the effective 75 for a doubled mode is equal to +75 or —75 
depending on the number of ir/a momentum components is even or odd. Since 
the sign of 75 is usually taken to be the sign of the chiral charge, then we have 
2 d ~ 1 — 1 doubled modes of positive chiral charge, but 2 d_1 of negative chiral 
charge. After taking into account of the primary mode, the number of fermion 
modes of positive chiral charge is equal to the number of negative ones. If all 
these modes contribute to the propagator, then the total chiral charge vanishes, 
the axial anomaly cancels, and the index is zero in any background gauge field. 
That is the case for the naive lattice fermion. 

For the Neuberger operator, the suppression of doubled modes depends on 
the factor T(p) which is a function of mo- At the origin and the corners of the 
BZ, sin^a) = 0, then T(p) can be simplified to the following, 

T{p) = \[\u{p)\-u{p)] (31) 

where the possible values of u(p) are the following 

u(p) = -TO 0) -m + 2r w , -to + 4r w , ■■■ , -m + 2dr w (32) 

for these 2 d momenta. From Eq. fl3~l|), it is evident that T(p) = for u(p) > 0. 

For m < 0, T(p) = for all 2 d massless fermion modes, hence they are 
all suppressed and decoupled from the theory. Since the theory does not have 
any massless fermion modes, so the index must be zero after the background 
gauge field is turned on. This is a topologically trivial phase of D. Note that 
for mo < 0, D can be shown to be topologically trivial for any background 
gauge field since D w does not possess any negative eigenvalue for a positive 
mass parameter. 

For all values of mo > 0, the primary mode at p = is recovered to the 
theory since u(0) < and T(0) ^ 0. For < m < 2r w , u(p) > and T{p) = 
for all doubled modes, so all doubled modes are still decoupled from the theory. 
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When we turn on the background gauge field with topological charge Q top , the 
primary mode then gives the index equal to +Qt op , same as the continuum 
theory. This is the topologically proper phase of D. 

At the phase boundary m = 2r w , all doubled modes have T(p) = 0, hence 
they are decoupled from the theory. However, when we turn on the background 
gauge field with topological charge Qt op , the index is always zero. From the 
viewpoint of perturbation theory, it is not clear what has happened. If we 
examine the doubled modes more closely, we find that there are d doubled 
modes which are just "critically" decoupled, i.e., each having only one nonzero 
( n /a ) momentum component which gives u{p) = and T(p) = 0, hence they 
are decoupled from the theory. On the other hand, if we assume that these d 
doubled modes still remain in the theory, then the index should be (1 — d)Q top 
rather than zero, since they carry negative chiral charge. So, this assumption 
cannot be true. Moreover, T{p) goes to zero smoothly at these d double modes. 
Hence, the theory is well behaved and free of massless doublers at the phase 
boundary m = 2r w in the free fermion limit but becomes topologically trivial 
in any background gauge field. This " anomalous" phenomenon seems to be of 
nonperturbative or topological origin since it cannot be understood from the 
viewpoint of perturbation theory. Further discussions of the behaviors around 
the phase boundaries are given in the Appendix. 

For 2r w < m < Ar w , there are d doubled modes each having only one 
nonzero ( n/a ) momentum component which gives u{p) < and T(p) ^ 0, 
hence they are restored to the theory and each one of them carries negative 
chiral charge. When we turn on the background gauge field with topological 
charge Qt op , these d double modes produce an index —dQ top which adds to the 
index +Qt op of the primary mode, thus gives the total index (1 — d)Q top . It 
agrees with numerical experiments for D in smooth background gauge fields 
except at the phase boundary, m = 4r w , where the actual index is always 
zero. 

For 4r w < m Q < 6r w ( provided that d > 4 ), there are additional d(d—l)/2 
doubled modes restored to the theory. Each of them has only two nonzero mo- 
mentum components which give u(p) < and T(p) ^ 0, hence carries positive 
chiral charge. Then the total chiral charge of the theory is 1 — d + d(d — l)/2 
after taking into account of the primary mode, the d doubled modes having 
one non-zero momentum component, and d(d — l)/2 doubled modes having 
two non-zero momentum components. When we turn on the background gauge 
field with topological charge Qtop, the index is equal to (2 — 3c? + d 2 )Q top /2. 
This is a topologically improper phase. 

The process of dividing into topologically improper phases each of width 
2r w continues until m = 2dr w , then for m > 2dr w , T{p) ^ for all doubled 
modes and they are all restored to the theory. The index must be zero for any 
gauge configurations, thus the theory is topologically trivial for m > 2dr w . 

Now we can pull the exact solutions of eigenvalues and the fermion prop- 
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agator together to obtain a complete picture of the topological phases in the 
free fermion limit. From the exact solutions of eigenvalues, Eq. ([24|), and the 
fermion propagator, Eqs. (EH), and (|5T]) , we immediately obtain that the 
+2 eigenmodes must have T(p) = in the fermion propagator, thus are com- 
pletely decoupled from the theory; while a zero mode has T(p) = \u(p)\, hence 
its presence depends on the value of mo. Therefore, using Eq. (p3|) , from the 
total chiral charge of the zero modes, we can obtain a formula for the index 
of D in a smooth backgroundgauge field. From Eqs. ([31])-([3^), we obtain the 
general formula for the total chiral charge of the zero modes in each phase, 

nW f { ~ iridl k-l...d(33) 

where the integer k denotes the k-th phase region, 2(k — l)r w < mo < 2kr w 
which includes the upper phase boundary, and Q$ = for m < and m > 
2dr w . Inserting (^) into Eq. (p3|), we obtain the index of D in each phase as 
the following 



index [D 



(*)l 



' if k = for m < 0, 

(-i) k+1 d ! Q if t - 1 ... W f or 
d r(i+d~fc) r(fe) ^*°p 11 h ' u iUi ^34) 

2(/c — l)r w < m < 2kr w , 
if k — d+ 1 for niQ > 2dr w . 



It is evident that (|34| ) satisfies the exact reflection symmetry index[I)( fc )] = 
— mdex[D^ d+1 ^ k ^] ( equivalent to Eq. ( |17D ) except at each one of the phase 
boundaries ( tuq = 2kr w , k — 1, • • • , d ) where the index is nonzero and equal to 
the index on the LHS of the boundary. Then it is obvious that (|34l) is invalid 
at the symmetry axis m = dr w where the index must be zero, Eq. (0), 
according to the exact reflection symmetry. In general, we can prove that fl34|) 
breaks down at each one of the phase boundaries ( m = 2kr w , k — 1, • • • , d ). 
Proof : 

Assuming Eq. (|34l) holds at the phase boundary m = x, then we have 

index[Z}(x)] = index[D(x — e)] (35) 

where e is an infinitesimal positive real number. Then, the image of tuq = x 
with respect to the mirror at mo = dr w is m = 2dr w — x, and its index is 

index [D{2dr w — x)} = — mdex[D(x)} (36) 



according to the exact reflection symmetry, Eq. (]17|) . Now consider the image 
of mo = x — e, and its index is 

index [D{2dr w — x + e)] = — index [D{x — e)] (37) 
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Then from Eqs. (|35|)-(|37|), we obtain 

index [D(2dr w -x+e)] = index [D(2dr w - x)] (38) 

This leads to the contradiction since the indices of two neighbouring phases 
must be different. Therefore Eq. ( |35| ) must be invalid since the exact reflection 
symmetry is always true under any circumstances. This completes the proof 
that (H]) fails at each one of the phase boundaries ( m = 2kr w , k = l,---,d 

)• 

In Table 1, we list the real eigenmodes of D in the free fermion limit on a two 
dimensional lattice with periodic boundary conditions. The lattice can be the 
infinite lattice or a finite lattice with even number of sites in each dimension. 
First, we recall that the zero modes and +2 eigenmodes of D are both chiral 
H, thus their chiral charges can be determined by the number of nonzero ( 
7r/a ) momentum components of this fermion mode. We also recall that for 
each real eigenvalue of D in two dimensions, it has degeneracy 2 which consists 
of a pair of eigenmodes of ±1 chiralities, while complex eigenvalues must come 
in complex conjugate pairs with chirality zero. The chirality degeneracy of 
real eigenmodes is not shown explicitly in Table 1. In the first phase, — oo < 
mo < 0, there is no zero modes and the theory in fact does not describe 
a massless fermion. The factor T(p) is zero at the origin p = as well as 
at the corners of the Brillouin zone. In the second phase, < mo < 2r w , 
there is one zero mode with positive chiral charge. The theory in this phase 
correctly describes a single flavor of massless fermion. At the phase boundary, 
mo = 2r w , two +2 real eigenmodes and their chirality partners become two 
pairs of complex conjugate eigenmodes ( 1 ± i ) which are decoupled from the 
theory, i.e., T(p) = 0. The number of zero modes remains the same. These two 
pairs of complex conjuagte eigenmodes become two zero modes ( in chirality 
pairs ) in the third phase, 2r w < m < 4r w . At the phase boundary, m = 4r„,, 
the number of zero modes remains the same, but the +2 eigenmode ( and 
its chirality partner ) become a complex conjugate pair ( 1 ± i ). Then they 
reappear as a zero mode ( in chirality pair ) in the last phase, 4r w < m < oo, in 
which all real eigenmodes are zero modes. We notice the interplay between the 
zero modes and the +2 eigenmodes when we change the value of m from — oo 
to +oo. Only in the region < mo < 2r w and the critical point mo = 2r w , the 
theory has the correct particle contents to describe a single flavor of massless 
fermion, while other regions suffer from either species doubling or absence of 
zero modes. The total chiral charge of zero modes in each phase is listed in the 
last column, from which we can predict the index of D in a smooth background 
gauge field according to Eq. (|23|), i.e., Eq. ([34]). As we will see, when we turn 
on a smooth background gauge field with non-zero topological charge, Eq. 
(PH) is satisfied exactly for all phases except at the phase boundaries where 



the index is always zero regardless of the total chiral charge ( see Table 3 ). 
In Table 2, we list the real eigenmodes of D in the free fermion limit on a 
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four dimensional lattice with periodic boundary conditions. Although it has 
many more phases than Table 1 for two dimensions, however, it shares the 
same essential features of Table 1, namely, only the region < m < 2r w and 
the critical point m = 2r w can describe a single flavor of massless fermion. 
We also recall that for each real eigenvalue of D in four dimensions, it has 
degeneracy 4 which consists of two pairs of eigenmodes of ±1 chiralities, while 
complex eigenvalues come in complex conjugate pairs with chirality zero and 
with degeneracy 2. The chirality degeneracy of real eigenmodes is not shown 
explicitly in Table 2. We will not go through each phase in Table 2 as what we 
have done for Table 1 since it is straightforward to obtain these results from 
the exact solutions of eigenvalues, Eq. and the fermion propagator, Eqs. 
©, © and ©. 

Our next task is to turn on the background gauge field to see to what extent 
the formula, Eq. (P4|), which relies on the weak coupling perturbation theory 
and the exact solutions in the free fermion limit, is satisfied in a background 
gauge field with nonzero topological charge, as well as to see the displacements 
of the phase boundaries with respect to the background gauge field. 



4 Topological phase diagrams 
for (1=2 and d=4 

In this section we turn on the background gauge field to investigate topo- 
logical phases of D on two dimensional and four dimensional lattices respec- 
tively. First, we define the background gauge fields in the following. On 



a 4-dimensional torus ( x M G [0, L^],^ = 1, 
SU (2) gauge fields can be represented as 

2nhi 2nq\x 2 



A 1 {x) 
A 2 (x) 
A 3 {x) 
AAx) 



+ 4 0) sin 



T 



2nh 2 
2irh 3 
27r/i 4 



I, the simplest nontrivial 
(2^n 2 



V u 



-x 2 



|(0) 



sin 



+ A 2 
2nq 2 Xi 



/27T77.1 



-Xi 



(0) 



3-^4 



+ A 



(0) 



sin 



+ A 



( 2Tm? 



sin 



/27m4 



-Xi 



V L, 



-x 3 



(39) 
(40) 
(41) 
(42) 



where r is one of the generators of SU(2) with the normalization tr(r 2 ) = 2, 
and qi and q 2 are integers. The global part is characterized by the topological 
charge 

Q = t^U / <?x e^xa tr(F^F A(T ) = 2q x q 2 (43) 
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Table 1: The real eigenmodes of D in the free fermion limit on a two dimen- 
sional lattice with periodic boundary conditions. The first column lists the 
ranges of m values, the second column the real eigenvalues, the third column 
the multiplicities of the eigenvalues in the momentum space, the fourth column 
the effective 75 of the eigenmode where the sign of 75 is used to denote the 
chiral charge of the eigenmode, and the last column is the total chiral charge 
of the zero modes. These results are obtained analytically from the exact so- 
lutions on a finite lattice of even number of sites in each dimension, or the 
infinite lattice. 
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Table 2: The real eigenmodes of D in the free fermion limit on a four dimen- 
sional lattice with periodic boundary conditions. The first column lists the 
ranges of m Q values, the second column the real eigenvalues, the third column 
the multiplicities of the eigenvalues in the momentum space, the fourth column 
the effective 75 of the eigenmode where the sign of 75 is used to denote the chi- 
ral charge of the eigenmode, and the last column the total chiral charge of the 
zero modes. These results are obtained analytically from the exact solutions 
on a finite lattice of even number of sites in each dimension, or the infinite 
lattice. 
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which must be an integer. The harmonic parts are parameterized by four 
real constants hi, h%, h 3 and h±. The local parts are chosen to be sinusoidal 
fluctuations with amplitudes Af\ A 2 °\ Af^ and A± , and frequencies 
^j^, ^jjf - and where n\, n 2 , n 3 and n 4 are integers. The discontinuity of 
Ai(x) ( A 3 (x) ) at x 2 = L 2 ( x 4 = L 4 ) due to the global part only amounts to 
a gauge transformation. The field tensors F\ 2 and F34 are continuous on the 
torus, while other F's are zero. To transcribe the continuum gauge fields to 
the lattice, we take the lattice sites at x M = 0, a, (iV M — l)a, where a is the 
lattice spacing and L M = N^a is the lattice size. Then the link variables are 



Ui(x) = exp[L4i(x)a] (44) 
U 2 (x) = exp 



iA 2 {x)a + id X2 ^ N2 -r) a — r 



U±(x) = exp 



iA 4 (x)a + iS X4 ^ N4 -i) a — - r 



(45) 
(46) 
(47) 



U 3 (x) = exp [L4 3 (a;)a] (46) 

2ivq2X 3 

L 3 

The last term in the exponent of U 2 (x) ( U±{x) ) is included to ensure that 
the field tensor F\ 2 ( -F34 ) which is defined by the ordered product of link 
variables around a plaquette is continuous on the torus. Similar construction 
can be carried out for the background gauge field on a two dimensional lattice, 
which has been described in details in [F3], and we use the same notations as 
defined in Eqs. (7)-(ll) in 0. 

In Table 3, we list the real eigenmodes of D on a 6 x 6 lattice in the 
background gauge field with topological charge Q = 1, and other parameters 
ht = 0.5, ha = 0.6, 4° } = 0.7, A { 2 0) = 0.8 and n x = n 2 = I. The Wilson 
parameter r w has been set to 1.0. The indices of D in the last column are 
computed by counting the zero modes according to index(D) = n_ — n + . They 
agree exactly with the formula, Eq. (|34]), with the total chiral charge of the 
zero modes obtained from the exact solutions in the free fermion limit as listed 
in the last column of Table 1, except at the phase boundaries where the index 
is always zero regardless of the total chiral charge. The response of D in the 
first region is easily understood since it has no zero modes in the free fermion 
limit and thus it cannot possess any zero modes after we turn on the gauge field 
with nonzero topological charge. The phase boundary has been shifted from 
m = to m = 0.5472 due to the background gauge fields. The index of D 
is always zero in this phase and D is called topologically trivial. In the second 
phase, the index of D is always equal to the topological charge of the gauge 
field, and D is called topologically proper. Since D is free of species doubling 
and the chiral charge of the zero modes is +1, thus D satisfies the Atiyah- 
Singer index theorem exactly in this phase. However, at the phase boundary 
of the central critical line mo = 2.0, D has the correct particle contents in 
the free fermion limit, but it does not possess any zero modes after turning on 
the topologically nontrivial background gauge field. This cannot be understood 
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real eigenvalues 
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Table 3: The index of D on a 6 x 6 lattice in the background gauge field with 
topological charge Q = 1, and other parameters hi = 0.5, h 2 = 0.6, = 0.7, 
4 0) = 0.8 and ni — n 2 = 1- The Wilson parameter r w is set to 1.0. We 
note that the gauge field shifts the critical points from their values in the free 
fermion limit. The reflection symmetry is satisfied exactly by checking that 
the indices are odd under the reflection with respect to the symmetry axis 
mo = 2.0, and the displacement at the first phase boundary 5q = 0.5472 is 
equal to minus of the displacement 5 2 = 3.4528 — 4.0 = —0.5472 at the third 
phase boundary, and the location of the symmetry axis m = 2.0 is invariant 
for any background gauge field. All critical points are accurate up to ±0.00005 
except the central critical point 2.0 is exact. The indices in the last column 
agree exactly with Eq. ( |3"4] ) except at the phase boundary ( m = 2.0 ). 



from the viwepoint of perturbation theory, though we know that the index must 
be zero at the symmetry axis, Eq. ([IB]), due to the exact reflection symmetry. 
Further discussions of the behaviors at the phase boundary are given in the 
Appendix. If we step into the next phase, we find that the index is opposite 
to the topological charge, in complete agreement with the formula (0) with 
the total chiral charge Q$ = — 1 in the free fermion limit, as listed in Table 
1. In the last phase, D has four species of massless fermion modes in the free 
fermion limit, similar to the case of the naive lattice fermion, then there is no 
zero modes and the index is zero for any background gauge field. We note that 
the displacements of the phase boundaries due to the background gauge field 
must obey the exact reflection symmetry 5k = —Sd-k, k = 0, ■ ■ ■ ,d, which is an 
immediate consequence of Eq. (|P7D. In Table 3, 8 = 0.5472 at the first phase 



boundary is equal to minus of the displacement 5 2 = 3.4572 — 4.0 = —0.5472 at 
the third phase boundary. The central phase boundary mo = 2.0 is invariant 
for any background gauge field. In Fig. 2, we sketch the topological phase 
diagram of D in smooth background gauge fields on a two dimensional lattice, 
which summarizes the gross features of a lot of numerical experiments. The 
discrete symmetry, index [D (4r w —m )] = — index [D(m )} is satisfied exactly in 
all cases. The index in each phase agrees exactly with the formula (|3"4]), with 
the total chiral charge of the zero modes obtained from the exact solutions in 



17 



A 



o 





/ 








/ 


| \ 






/ 

/ 


1 \ 

1 \ 
1 \ 






/ 


1 \ 

1 \ 




trivial 


/ proper 


| improper \ 


trivial 






1 \ 




index = 


/ index = Q 


| index = -Q \ 


index = 




/ 


1 \ 









2 r. 



4 r. 



Figure 2: The topological phases of D on a two dimensional lattice for smooth 
background gauge fields. The phase boundaries denoted by dashed lines are 
sketched to indicate their evolutions with respect to the background gauge 
field. The reflection symmetry with respect to the central critical line mo = 2r w 
is displayed explicitly. Note that the central critical line m = 2r w is invariant 
for any background gauge configurations. The phase transition points in the 
free fermion limit ( at the horizontal line ) are determined exactly. 



the free fermion limit as listed in the last column of Table 1, except at the 
phase boundaries where the index is always zero regardless of the total chiral 
charge. 

In Table 4, we list the real eigenmodes of D on a 4 4 lattice in the background 
gauge field with topological charge Q = 2, and other parameters hi = 0.1, 
h 2 = 0.2, h 3 = 0.3, h A = 0.4, A\ 0) = 4 0) = 4 0) = 4 0) = 0. The Wilson 
parameter r w is set to 1.0. There are more phases regions than those in two 
dimensions ( see Table 3 ), and the locations of all critical points are shifted 
by the background gauge field except the central critical point m = 4.0. It 
is evident that the reflection symmetry, Eq. (|ITD is satisfied exactly. We note 
that the displacements of the phase boundaries due to the background gauge 
field are 5 = 0.3762, 8 X = 2.1881 - 2.0 = 0.1881, 5 2 = 0, 5 3 = 5.8119 - 
6.0 = -0.1881 and 5 4 = 7.6238 - 8.0 = -0.3762, which satisfy the exact 
reflection symmetry 5k = —5^^,^ — 0, •■•,4. We do not go through the 
detailed descriptions for each phase as we have done for Table 3. We have 
also performed extensive numerical tests in the topologically proper phase as 
well as in the improper phases. For smooth background gauge fields with 
nonzero topological charge, the exact zero modes with definite chirality are 
reproduced to a very high precision on a finite four dimensional lattice, and 
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Table 4: The index of D on a 4 4 lattice in the background gauge field with 
topological charge Q = 2, and other parameters hi = 0.1, hi = 0.2, h 3 = 0.3, 
h 4 = 0.4, AP = A { 2 ] = = Af } = 0. The reflection symmetry, index[D(8- 
mo)] = — index [D(m )} is satisfied exactly, in particular, the displacements of 
the phase boundaries due to the background gauge field satisfy the relationship 
5k = —S^kyk = 0,---,4. All critical points are accurate up to ±0.00005 
except the central critical point 4.0 is exact. The indices in the last column 
agree exactly with the formula (|34|), except at the phase boundaries m = 
2.1881, 4.0, 5.8119. 
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Figure 3: The topological phases of D on a four dimensional lattice for smooth 
background gauge fields. The phase boundaries with dashed lines are sketched 
to indicate their evolutions with respect to the background gauge field. The 
reflection symmetry with respect to the central critical line mo = 4r w is shown 
explicitly. Note that the central critical line mo = ^r w is invariant for any 
background gauge configurations. The critical points in the free fermion limit 
( at the horizontal line ) are determined exactly. 



the index theorem is satisfied exactly in the topologically proper phase. The 
zero modes are also very stable under random fluctuations as defined in Eq. 
(43) of ref. 0. This completes our intended investigations for D in four 
dimensional smooth background gauge fields as outlined in ||. In Fig. 3, we 
sketch the topological phase diagram of D in smooth background gauge fields 
on a four dimensional lattice with even number of sites in each dimension, 
which summarizes the gross features of a lot of numerical experiments. The 
discrete symmetry, index[D(8r w — mo)] = — index[D(mo)] is satisfied exactly 
in all cases. The index in each phase agrees exactly with the formula, Eq. 
(|34"D, except at the phase boundaries where the index is always zero regardless 
of the total chiral charge. 

It is interesting to note that in ref. |J, the spectral flow of the her- 
mitian Wilson-Dirac operator H w (mo) = / j 5 D w (m ) is studied numerically 
for smooth SU(2) instanton backgrounds on 8 4 lattice, and the symmetry 
—H w (8 — m , U) = H w (m , —U) is displayed. Since the index of D can also 
be computed by studying the level crossings in the spectral flow of H w as 
a function of m [ see Eq. (^1| ) in the Appendix ], the results of ref. 
provide another four- dimensional verification of the reflection symmetry (17) 
proved in Section 2. Furthermore, in ref. [110], the spectrum of the hermitian 
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Table 5: The index of D on a 6 x 6 lattice in the background gauge field with 
topological charge Q — 1, and other parameters h\ = 0.5, h 2 = 0.6, = 2.0, 
= 0.8 and ri\ = n 2 = 1. The Wilson parameter r w is set to 1.0. The 
topological phases have bifurcated from 4 into 13 regions in two steps. The 
reflection symmetry with respect to the central critical line tuq = 2.0 is satisfied 
exactly by the indices in all phases as well as the positions of the critical points. 
All critical points are accurate up to ±0.00005 except the central critical point 
2.0 is exact. 



Neuberger-Dirac operator H = 75D in a smooth SU(2) instanton background 
is investigated in the phase < itlq < 2. The chiral properties of 77 ( 77 ) 
are verified and the pole method 



11 1 appears to be the most cost effective 



implementation. 



5 Summary and Discussions 

We have investigated the topological phase diagrams of the Neuberger-Dirac 
operator D with respect to the negative mass parameter mo and the back- 
ground gauge field. When the gauge configurations are far from being rough, 
the phase diagrams can be schematically represented in Figs. [I], ||] and [3] for 
(i-dimensions, 2-dimensions and 4-dimensions respectively. Each topological 
phase diagram ( Figs. 1-3 ) possesses the exact reflection symmetry (|T7|). The 
symmetry axis mo = dr w is a phase boundary with zero index and its posi- 
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tion is invariant for any background gauge field. Among the many different 
phases in each phase diagram, there is only one phase in which the Atiyah- 
Singer index theorem is satisfied exactly, while other phases are either trivial 
or improper. 

The formula (Q) for the index of D in each phase is derived by obtaining 
the total chiral charge of the zero modes in the exact solution of the free fermion 
propagator. For smooth gauge configurations, (|34f) is verified to be correct by 
numerical experiments on two and four dimensional finite lattices with even 
number of sites in each direction, except at the phase boundaries starting at ( 
mo = 2kr w , k — 1, • • - , d ) where the index is always zero for any background 
gauge fields. Moreover, using the exact reflection symmetry, we prove that 



( |34"D must break down at the phase boundaries. This indicates that there 
may be some nonperturbative or topological effects happening at the phase 
boundaries. Further discussions of the behaviors at the phase boundaries are 
given in the Appendix. Now it is natural to generalize ([34] ) by incorporating 
the displacements of the phase boundaries due to the background gauge field, 
and also imposing the index zero at all phase boundaries, then fl34j) becomes 

{ d r(i+d-k) r(fe) Qt°p if 2 ( k ~ l ) r ' w + <m o< 2kr w + 5 k 
for k = 1,- --,d; 
if m = 2kr w + 5^ for k — 0, • • • , d. 

(48) 

where 6k, k = 0, • ■ ■ , d are the displacements of the phase boundaries due to 
the background gauge field, which satisfy 

4 = S d - k . (49) 

Equations fl48|) and ([|9|) constitute the general formula for the index of Neu- 
berger operator in a smooth background gauge field, as a function of the neg- 
ative mass parameter m^. The exact reflection symmetry ( |i~7l) is obviously 
satisfied by Eqs. and (gDp. 



It is tempting to go beyond the limit of smooth background gauge field 
by increasing the field strength or its roughness, in order to explore what lies 
ahead in the upper portions of these phase diagrams. However, we quickly run 
into difficulties as the phase boundaries bifurcate into many different phases, 
reminiscent the chaos in dynamical systems. Even for D in a topologically 
proper phase, an infinitesimal change in mo may end up having D in a trivial 
or improper phase. For a given rough gauge configuration, we do not know 
a priori what is the range of mo values for D to be in the topologically 
proper phase. We can easily demonstrate the bifurcation of the topologcial 
phases of D by increasing the amplitudes of the local sinusoidal fluctuations 
of the background gauge field on a two dimensional lattice. In Table 5, all 
parameters are the same as Table 3 except the amplitude is increased 
from 0.7 to 2.0. We see that there are total 13 phase regions which are resulted 
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from bifurcations in two steps. It is evident that formula (pESD fails completely. 
However, the exact reflection symmetry, Eq. fllTD , for the index in each phase 
as well as the positions of the phase boundaries, is still satisfied exactly, though 
we cannot predict the number of different topological phases and the index of 
D in each phase. We do not know precisely the conditions under which the 
bifurcations may happen, but one thing for sure is that the index is identically 
zero at all phase boundaries. 

We have also measured the indices of D with mo fixed at 1.0 for an ensemble 
of gauge configurations from quenched simulations, and compared them to 
the geometrically defined topological charges of the gauge configurations. It 
turns out that the index does not always agree with the geometrically defined 
topological charge. The discrepancies may be due to the fact that the fixed 
value m = 1.0 does not fall into the topologically proper phase for all gauge 
configurations. If that is the case, we can proceed in the following way : for 
a given gauge configuration, we measure the index of D(mo) for sufficiently 
sampled itlq values and exploit the reflection symmetry, then we should be able 
to identify the correct index unambiguously. Otherwise, we must resolve the 
discrepancies from the fundamental viewpoint, namely, the very definition of 
topological charge on a lattice. 

Finally we would like to remark that it is interesting to understand the 
underlying dynamics which provoke the collective motion of the eigenvalues of 
D ( or equivalently the zero crossings of the spectral flow of the hermitian 
Wilson-Dirac operator H w = ) at the phase boundaries ( except for the 

first and the last phase boundaries ) such that the index always turns out to 
be zero, though the exact refection symmetry must impose the zero index at 
the phase boundary m = dr w . 



In this appendix, we discuss the spectral flow of the Neuberger operator as a 
function of the negative mass parameter m , in particular, its behaviors at the 
phase boundaries, in order to gain a perspective on the underlying mechanism 
which gives the zero index at the phase boundaries. Note that the index at 
the phase boundary m = dr w ( the symmetry axis ) must be zero due to the 
exact reflection symmetry (|P71). Since the Neuberger operator can be written 
in terms of the hermitian Wilson-Dirac operator H w (m ) = r j^D w {mo) ) 

D = 1 + V = 1 + 75-7==, (50) 
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then the index of D is 

index (D) = ±Tr( l5 V) = ^Tr ^k=j = \{h + - h.) (51) 

where Tr denotes the trace over the Dirac space and the position space, h + {h ) 
is the number of positive ( negative ) eigenvalues of H w . So, the index of D 
can be obtained by studying the level crossings of the low-lying eigenvalues 
of H w M. However, the spectral flow of the Neuberger operator D seems to 
be more complicated than that of the hermitian Wilson-Dirac operator H w . 
Since eigenvalues of H w (m ) are real, its spectral flow as a function of m 
can be plotted as curves on a plane. On the other hand, the eigenvalues of 
the Neuberger operator are complex, lying in the unit circle with center at 
one on the complex plane ||, therefore the spectral flow as a function of mo 
can be plotted as curves on the unit cylinder with center at one in the three 
dimensional space formed by the complex plane of the eigenvalues and the 
parameter mo- The relationship between the eigenvalues of D and H w can be 
written as 

A = 1 + U*j 5 U H ^Ltf H U (52) 

/\2 
V A H 

where A ( A# ) is the diagonal matrix of eigenvalues of D ( H w ), U ( Uh ) is 
the unitary transformation matrix formed by the eigenfunctions, satisfying 

D = U\U ] (53) 
H w = U H X H Ui (54) 

It is evident that the spectral flow of D not only depends on the spectral flow 
of H w , but also the overlap of their eigenfunctions. 

In general, due to the limited precisions of most numercial computations, 
one cannot trace the spectral flow of D to any precisions one wishes. However, 
even if one can trace the spectral flow to a very high precision, one still do 
not understand why the eigenvalues behave in such a manner at the phase 
boundaries so that the index must be zero, except for the first and the last 
phase boundaries. For example, let us describe the spectral flow of D around 
the phase boundary tuq = 2 in Table 3. Although we know that its index 
must be zero due to the exact reflection symmetry, it is still instructive to 
examine the spectral flow around this phase boundary. At rriQ = 2 — e, where 
e is any positive infinitesimal number, the zero mode of negative chirality 
is accompanied by a +2 eigenmode of positive chirality, and the rest of the 
eigenvalues are in complex conjugate pairs. However, at m = 2, the zero 
mode becomes a member of a complex conjugate pair ( a\ ± ia 2 ), while the 
associated +2 eigenmode becomes a member of another complex conjugate pair 
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( 2 — cz3 ± id4 ), where aj(i = 1, • • • , 4) are positive small numbers. Presumbaly 
these two pairs of complex conjugates have no direct correlations since they are 
far apart. However, in order to create one additional complex conjugate pair 
out of the zero mode and its associated +2 real eigenmode, all other complex 
conjugate pairs must undergo a collective motion in the unit circle such that 
all eigenvalues turn out to be in complex conjugate pairs. Then at m = 2 + e, 
all complex conjugate pairs again undergo another collective motion to yield 
a zero mode of positive chirality and its associated +2 eigenmode of negative 
chirality. It seems to be unlikely that such collective motions of the eigenvalues 
at the phase boundaries can be described by any perturbation theory in itlq or 
the background gauge fields. 

Now let us turn to the spectral flow of the hermitian Wilson-Dirac operator 
H w = 75-D™, which corresponds to the spectral flow of D described above. At 
m = 2 — e, the eigenvalues of H w are mostly in plus/minus pairs except two 
positive small eigenvalues. Thus h + — h_ = 2, and the index of D is one 
according to Eq. (|5TD. At m = 2, one of the positive small eigenvalues crosses 
zeroQ and becomes a negative small eigenvalue. The transition is exactly at 
the point m = 2. So h + = /i„, and the index is zero. Then at m = 2 + e, the 
remaining positive small eigenvalue also crosses zero and turns into a negative 
small eigenvalue. Since h_ — h + = 2, the index is —1. Now the seemingly puzzle 
is why these two positive small eigenvalues do not simultaneously become two 
negative small eigenvalues at the point m = 2, then the phase boundary of 
zero index would not exist and the formula ([34]) would also hold for the phase 
boundaries. However, this would lead to a contradiction to the exact reflection 
symmetry (|18|) which must hold in general, and the only way out is to have the 
spectral flow to produce zero index at the phase boundary of the symmetry 
axis ( m = dr w ). So, this puzzle seems to be resolved, but the underlying 
mechanism is still not understood. In general, for phase boundaries other than 
the first one ( stemming at m = ) , the symmetry axis mo = dr w , and the last 
one ( stemming at m = 2dr w ), we do not have any arguments to explain why 
there must be the number ^\h+ — h-\ of eigenvalues undergo a sign-flipping 
transition when approaching one side of the phase boundary, while the number 
h\h' + — h'_\ the other side of the phase boundary, such that the index must be 
zero at the phase boundary. Furthermore, it seems to be more difficult to 
tackle what is the underlying dynamics that provoke such transitions at the 
phase boundaries. 



2 If one examines the crossing numerically, one may find a seemingly discontinuity of order 
10 , however, which is only due to the limited precision of one's numerical computation. 
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